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NOMENCLATURE 


F ,  G  ,  H 
i  ,  j  ,k 
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T  ,  T  ,  T 
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time 
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mass  averaged  dependent  variable 
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density 
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SECTION  1 
INTRODUCTION 


The  assessment  of  unsteady  flow  structure  around  an  airborne 

1  2 

turret  is  of  current  interest  in  the  aero-optical  research  area  '  . 

The  reduction  of  cumulative  total  system  optical  error  is  becoming 
increasingly  important  for  design  purposes.  The  major  contributor 
to  such  degradation  is  the  fluid  mechanical  interference.  In 
particular,  the  complex  turbulent  mixing,  entrainment  process  and 
the  separated  flow  result  in  free  shear  layer  oscillations .  Although 
a  complete  physical  understanding  of  this  fluid  dynamics  is  lacking, 
significant  progress  has  been  made  in  the  past  few  years. 

To  estimate  the  aerodynamic  induced  density  variation,  the 

research  efforts  were  directed  towards  using  experimental,  analyti- 

3  4 

cal  and  numerical  methods.  The  experimental  studies  '  addressed 
the  use  of  extensive  instrumentation  to  statistically  quantify  the 
turbulence  levels  and  to  provide  flow  visualization.  These  experi¬ 
mental  parametric  studies  were  conducted  for  various  flow  conditions. 
The  primary  objective  was  to  estimate  the  changes  in  the  flowfield 
fluctuations  as  a  function  of  small  scale  turbulence  in  the  shear 
layer  at  a  selected  station.  Meanwhile,  the  experimental  scaling 
laws’’  were  formulated  for  RMS  density  fluctuations,  transverse 
scale  length  and  shear  layer  width  with  state-of-the-art  unsteady 
aerodynamics  and  optical  measurements.  The  data  for  the  full  scale 
model  was  simulated  using  a  compressible  boundary  layer  program. 
However,  the  numerical  modeling  of  such  three-dimensional  problems 
is  prefered  because  it  provides  a  more  detailed  description  of  the 
flowfield  at  a  lower  cost.  In  order  to  arrive  at  the  fruitful 
results,  the  quantitative  analysis  of  wake  structure,  separated 
flow  and  its  domain  of  influence  demand  a  systematic  approach. 

Our  earlier  attempt^  was  aimed  at  validating  the  computational 
procedure  and  reconfirming  some  of  the  fundamental  results  for 
flow  past  a  large  protuberance  (turret) .  We  were  successful  in 
demonstrating  the  usefulness  of  a  hypothesis  that  the  skin  friction 
lines  on  the  surface  are  said  to  behave  as  a  continuous  vector 


1 


field  and  the  singular  points  in  this  field  obey  the  topological 
7  8 

rule  '  .  This  verification  also  confirmed  the  primary  separation 

line  and  horseshoe  vortex  structure  for  the  surface  mounted 

9 

obstacle,  boundary  layer  interaction  .  Complementing  the  efforts, 
the  purpose  of  the  present  investigation  is  to  analyze  the  wake 
region  and  to  estimate  the  effect  of  flow  control  on  the  separated 
flow  and  the  shear  layer.  The  computation  is  performed  using  the 
CRAY-1  computer  for  a  free  stream  Mach  number  0.6  and  Reynolds 
number  12.8  x  106/m.  The  compressible  Navier-Stokes  equations  are 
solved  for  flow  around  a  hemispherically  capped  turret  mounted  on 
a  flat  surface.  The  turret  has  cutout  at  an  azimuth  angle  120° 
from  the  free  stream  direction.  The  analysis  is  done  for  two 
cases:  when  the  turret  is  treated  as  (1)  a  solid  shell  and  (2)  a 

porous  shell  with  uniform  suction  rate  of  0.136  kg/s.  The  experi¬ 
mental  conditions  for  the  mass  suction  rate  were  duplicated  by 
imposing  a  negative  normal  velocity  component  with  a  magnitude 
of  3.7  m/s  over  the  porous  turret  surface.  The  aim  is  to  evaluate 
the  basic  fluid  mechanical  difference  in  the  flowfield  due  to 

3  4 

mass  suction  and  compare  the  results  with  the  experimental  data  ' 


SECTION  2 

GOVERNING  EQUATIONS 

For  classical  fluid  mechanical  problems,  the  Navier-Stokes 
equations  constitute  a  very  accurate  mathematical  model.  With 
adequate  inclusion  of  turbulence  terms,  it  provides  the  basic 
understanding  for  the  viscous-inviscid  interaction  and  separated 
flow  at  high  Reynolds  number.  However,  due  to  computing  speed 
and  memory  limitation,  certain  approximations  are  still  needed 
to  use  these  equations1^. 

Applying  the  fundamental  principles  of  mechanics  and  thermo¬ 
dynamics  for  a  continuous  media  in  a  Galilean  (absolute)  frame 
of  reference,  the  integral  forms  of  the  conservation  law  written 
in  the  following  form'*'*': 

Conservation  of  mass : 


/  p‘ 


dv  =  0 


Conservation  of  momentum: 


J"  p*U*dv  =  j'  o*n*ds+  J  p*f? 


Conservation  of  energy: 


J  p*E*dv  =  J'  (u*a-q)  *n*ds  +  J  p*fe*U*dv 


where  n  is  the  unit  normal  vector  on  the  surface  9v^.  bounding 
the  volume  vfc.  These  integral  forms  are  supplemented  with  the 
constitutive  relations  for  the  stress  tensor  if  and  the  heat  flux 
vector  q  in  terms  of  the  basic  variables: 


(4) 


0  =  -p •  I  +  T 

q  =  -k • grad  T  ( f 

-* 

where  I  is  the  unit  tensor.  In  these  relations,  it  is  implicitly 
assumed  that  the  variables  p,  U  and  e  are  continuous  and  differ¬ 
entiable  . 

For  turbulent  flows,  which  characterize  the  large  variation 
in  fluid  physics,  the  numerical  solutions  of  the  equations  (1) 
through  (3)  is  feasible  only  if  turbulence  properties  are  replaced 
by  some  algebraic  relations  (turbulence  model)  and  the  equations 
are  written  in  a  more  suitable  form.  Thus,  in  the  absence  of  body 
forces,  the  unsteady,  three-dimensional  Navier-Stokes  equations, 
in  mass  averaged  variables  and  in  surface  oriented  coordinate 
system,  can  be  written  in  the  convenient  form  as  follows: 


3t  +  35  +  Cnx'ny,r|z] 


+  Px'S'^zJ 


where 


pu2+o 

X 

pUV+Txy 

PUW+Txz 

(Pe+cx)u+Txy.v+Txz*w-qx 


-+■ 

G  = 


r 

PV 

pVU+Tyx 

pv2+ay 

pvw+Tyz 

(pe+ay)v+Tyx-u+TyZ*w-4y 


pw 

pwu+t2x 

P"V+Tz 

pW2+Oz 

(pe+oz)w+Tzx.u+Tzy-v-qz 


(9) 


(10) 


X 


p  +  •-  (u+e) 


/ 3u  ,  3v  ,  3w\ 

\dZ  +  3^  +  Jz) 


2(P+c)-|S 


(ID 


Txy  =  V  “  '  <U+C)  (I?  +  ti) 

'xx  "  Tzx  "  ~  <u+£)  (If  +  If) 


^2  .  .  .  /3u  .  3v  ,  3w\ 
ay  =  P+3  (p+e)  +  ^  + 


-2(p+e).g 


(14) 


,  ,  .  /3v.  3w\ 

Tyz  Tzy  ~  y  £  \^z  3y/ 

2  ,  ,  ,  /3u  3v  ,  3w \  ,  3w 

cz  =  P  +  3  (y+e)  \33?  +  W  '2(y+E)  3? 


<*z  = 


Y(i?r  + 

y(pF  + 

Y(p?  + 


3e . 

_ i 

3z 


(18) 


(19) 


(20) 
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(21) 


o  =  +  (u2 +v2 +w2 ) /2 . 0 

e.  =  C  *T  (22) 

IV 

For  these  equations,  we  assume  that  the  thermodynamic  properties 
of  the  fluid  are  governed  by  the  perfect  gas  law, 

p  =  pRT  (23) 

with  constant  specific  heat  ratio,  y  =  1.4  (air)  and  the  coefficient 
of  viscosity  varies  with  temperature  according  to  the  Sutherland's 
viscosity  formula.  It  is  also  implied  that  the  thermal  conductivity 
coefficient  <  and  the  viscosity  coefficients  \,y  are  functions  of 
temperature  only  and  A,p  comply  with  the  Stoke 's  relation.  The 
molecular  and  turbulent  Prandtl  numbers  are  assumed  to  be  constant 
and  are  equal  to  0.72  and  0.90,  respectively. 

The  system  of  equations  (6)  through  (23)  is  closed  by 
assigning  a  specific  relation  for  the  coefficient  of  eddy  viscosity. 

The  fluid  motion  described  by  the  aforementioned  nonlinear, 
strongly  coupled  equations  is  to  be  analysed  around  a  hemispher- 
ically  capped  turret  (Figure  1)  mounted  on  a  flat  surface.  Its 
base  is  a  circular  cylinder  with  diameter  0.127  meter  and  height 
0.078  meter.  The  hemispherical  dome  at  the  top  has  a  diameter  of 
0.127  meter  whose  center  is  on  the  turret  axis.  If  we  measure 
the  (azimuth)  angle  clockwise  from  the  oncoming  stream  direction, 
the  turret  has  a  cutout  at  120°  AZA  at  a  normal  distance  0.057  meter 
from  the  turret  axis.  The  origin  of  the  coordinate  system  is 
located  at  the  intersection  of  the  turret  axis  and  the  flat 
surface.  The  £,n ,  S  directions  correspond  to  axial,  normal  and 
tangential  transformed  coordinates  of  the  turret,  respectively. 

The  y-direction  corresponds  to  the  turret  axis  in  the  Cartesian 
coordinates  (x,y,z) . 

One  of  the  controlling  factors  for  the  numerical  simulation 
is  the  proper  selection  of  coordinates  and  grid  system.  In  order 
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to  avoid  any  interpolation  in  the  transformed  computational  domain, 

the  effective  way  is  to  opt  for  a  turret  surface  oriented  grid 

system.  Since  such  a  system  has  coordinate  lines  coinciding  with 

the  surface  contours  of  the  turret,  all  boundary  conditions  can 

be  expressed  at  the  grid  points.  The  normal  derivatives  at  the 

turret  can  be  represented  using  only  finite  difference  between  grid 
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points  on  the  coordinate  lines  .  Thus,  we  define  the  turret 
surface  by  62  x  30  points  (Figure  2) .  On  the  cylindrical  portion, 
15  points  were  used  in  the  axial  direction  stretched  exponentially 
from  the  base  and  15  points  were  uniformly  placed  on  the  hemi¬ 
spherical  dome.  In  the  tangential  direction,  62  equidistant  points 
were  used  with  first  and  last  two  overlapping  points  to  permit 
the  use  of  fourth  order  pressure  damping.  Using  these  surface 
nodes  as  the  reference  points,  the  normal  coordinates  are  then 
described  by  30  exponentially  stretched  field  points  (x,y,z)  from 
the  turret  surface  extending  outwards  to  seven  times  the  turret 

diameter  (Figure  3) .  These  points  are  generated  by  the  homotropy 
13 

scheme  : 


(x,y,z) 


(x,y,z) 


outer 


(x,y,z) 


turret 


(24) 


where  C„  is  a  nonzero  constant.  The  finer  mesh  near  the  turret 
and  flat  plate  help  resolve  the  viscous  effects  while  the  coarser 
mesh  helps  to  reduce  the  computational  time. 

Once  the  coordinates  are  generated,  the  derivatives  of 
coordinate  transformation  can  be  obtained  through  the  Jacobian,  J 
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Figure  2.  Mesh  System  of  the  Turret 
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and  its  inverse 


j-1  _  3 (x,y , z)  _ 
3  ( S » H  /  4 ) 
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Care  has  been  taken  that  no  mesh  point  lies  on  the  turret  axis. 

This  precaution  allows  the  existence  of  the  inverse  Jacobian 
and  permits  a  one-to-one  and  onto  correspondence  from  the  physical 
to  the  transformed  space. 

For  this  configuration  and  the  mesh  system,  our  aim  is  to 

treat  the  turret  as  a  solid  shell  and  a  porous  shell  (with  uniform 

mass  suction) .  The  basic  input  for  the  distributed  suction  porous 
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shell  model  is  based  on  the  data  provided  by  Haslund  and  Reyhner 
The  maximum  flow  rate  available  from  the  suction  pump  of  0.136  kg/s 
was  shared  as  0.078  kg/s  through  the  hemisphere  portion  and 
0.058  kg/s  through  the  upper  half  of  the  cylinder  portion  (excluding 
cutout).  At  a  Reynolds  number  of  12.8  *  106/m  and  a  free  stream 
Mach  number  0.6,  this  suction  rate  eliminates  recirculation  between 
holes15.  For  our  mesh  system,  we  assume  that  the  suction  is  being 
applied  to  all  the  grid  point  locations  except  at  the  lower  half 
of  the  turret  and  the  cutout  area.  The  suction  Mach  number  cor¬ 
responding  to  the  assigned  suction  rate  is  0.01  which  results  in 
a  negligible  change  in  turret  surface  density  values. 

The  well  posed  and  stable  boundary  conditions  for  subsonic 
Navier-Stokes  equations  are  still  uncertain.  The  primary  physical 
phenomenon  consisting  of  the  propagation  of  pressure  waves  and 
vortices  limit  the  flow  variables  to  the  asymptotic  values  in  the 
farfield.  However,  for  three-dimensional  flows,  the  relief  effect 
is  significant  and  the  disturbances  generated  by  the  turret  decay 
rapidly.  Thus,  if  the  boundary  conditions  are  specified  at  a 
sufficiently  large  distance  from  the  turret,  the  overspecification 
is  generally  acceptable. 


1 1 


Initially,  we  assume  that  the  turret  is  immersed  completely 
in  the  uniform  stream.  The  upper  boundary  conditions  are  the 
free  stream  values,  whereas  at  the  downstream  side,  the  gradient 
of  the  flow  properties  is  assumed  to  vanish.  On  the  solid 
surfaces,  the  noslip  condition  for  the  velocity  component  along 
with  prescribed  wall  temperatures  are  imposed,  while  near  the 
turret  axis,  an  interpolation  is  adopted  to  eliminate  calculation 
on  this  singularity.  Along  the  coinciding  plane  (azimuth  angles 
0°  and  360°),  the  conditions  were  matched  at  all  the  mesh  points. 
Or,  in  other  words,  the  following  initial  and  boundary  conditions 
were  applied  on  the  unit  cube  (Figure  4)  in  the  transformed 
space  0  £  £,n,£  £  1.0. 

Initial  Conditions 


U(0,£,n,C)  = 


(27) 


Boundary  Conditions 


Upstream:  U(t,£,l,£0)  =  UM, 


(28) 


where  is  a  particular  value  of  C  corresponding  to  all  the  grid 
points  facing  the  upcoming  flow  (^  corresponds  to  the  azimuth 
angle  90°  and  for  the  azimuth  angle  270°). 


Downstream: 


3n 


=  0, 


n=l 


(29) 


On  the  solid  surfaces:  u  =  v  =  w  =  0 

T  =  1 .0695*T 
w  » 

3n 


(30) 


where  n  is  the  direction  normal  to  the  solid  surfaces. 


Figure  4.  The  Boundary  Conditions 
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(31) 
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Near  the  turret  axis:  U(t,l,n,c)  =  IPF 
where 

ds_  dSp+2  ds. 

IPF  =  - - -  .  u(t,l-AC,n,  C+30*AO  +  - — 

2(ds^+ds2)  2(ds^+ds2) 

x  U(t,l-A^,n,C) .  (32) 

It  may  be  recalled  that  there  are  60  equidistant  mesh  points 
in  the  tangential  direction  and  C  +  30*Ac  will  identify  the  diagonally 
opposite  mesh  point  across  the  turret  axis.  The  terms  ds^  and  ds2 
are  the  arc  lengths  subtended  by  the  mesh  points  (l,n,C)  and 
(l-A£,n»C)  along  the  turret  axis  to  the  center  of  the  hemisphere 
(turret  top) . 

These  are  the  boundary  conditions  when  the  turret  is  treated 
as  a  solid  shell.  For  the  computation  with  mass  suction,  the 
boundary  conditions  at  the  porous  turret  surface  are  replaced  by 
the  negative  normal  velocity  component  of  magnitude  3.7  m/s. 
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SECTION  3 
TURBULENCE  MODEL 


The  algebraic  eddy  viscosity  model  used  in  this  study  is  due 
to  Baldwin  and  Lomax^.  The  model  is  tailored  particularly  to 
suit  the  wake  region  at  the  turret  shoulder  where  substantial 
variation  in  all  the  three  directions  dominates.  It  is  described 
briefly  as  follows: 

In  the  inner  region,  the  eddy  viscosity  is  given  by 
^  =  (k1-L*D1)  2-p*  |oi|  (33) 

where  k^  =  0.4  is  the  von  Karman's  constant,  is  the  van  Driest' s 
damping  factor. 


and  oj  is  the  vorticity  function  (Appendix)  .  The  scale  length 
L,  is  the  distance  normal  to  the  turret  or  the  flat  plate  surface 
(Figure  5)  : 


L  =  min 


f(x-xT)2  +  (y-yT)  2  +  (z-zT)2 
<y-yFL) 


where  (x,y,z)  are  the  Cartesian  coordinates  of  the  field  point, 

(x,y,z)T  and  (x,y,z)FL  are  the  projection  coordinates  on  the  turret 

and  flat  plate  surface,  respectively.  This  relation  agrees  with 
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the  asymptotic  length  scale  formulation 
In  the  wake  region: 

eQ  =  Cc*Ca*FWAKE*FKIF.  (36) 

The  Clauser  contant  Cc= 0.0168  is  multiplied  by  an  additional 
constant,  Cft = 1 . 6  and  the  coefficient,  FWAKE,  is  given  by 
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f 

I 


FWAKE 


min 


L  •  TF 
max  max 


0 . 25  *L 


max 


ud2/tf 


max 


(37) 


The  scale  length  L  is  a  particular  value  of  L  when  the  function 

max 

TF  (=L*  I  co  I  )  attains  its  maximum,  TF  .  The  introduction  of  the 
11  max 

function  TF  in  this  turbulence  model  poses  certain  difficulties 

in  uniquely  defining  the  maximum  value  or  "the  peak".  In  the 

turret  wake  region,  this  function  strongly  depends  on  vorticity. 

Due  to  large-scale  oscillatory  flow,  the  function  TF  is  expected 
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to  have  more  than  one  peak  '  .  These  peaks  display  a  "double¬ 

hump”  like  character  causing  an  abrupt  change  in  L  .  Thus, 
during  our  computation,  we  always  take  the  first  peak  to  avoid  an 
inconsistency  and  possible  mutilation  in  the  final  results.  The 
other  terms  included  are  the  Klebanoff  intermittancy  correction 
factor,  FKIF . 


FKIF 


1.0  +  5.5 


_ L_ 

Lmax 


-1 


and  the  total  velocity  differential,  UD 


(38) 


UD  =  ^u2+  v2+ 


w 


-  ^u2+  v2+  w2 


max 


min 


The  effective  viscosity  is,  then,  taken  to  be 


e  =  min (e^ , eQ) . 


(39) 


(40) 
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SECTION  4 

NUMERICAL  PROCEDURE 


Owing  to  the  large  memory  availability  and  high  data  processing 
rate  (the  CPU  time  per  grid  point  per  time  iteration) ,  the  computer 
CRAY-1  was  chosen  for  our  computation.  The  well-established  vector 
solver^  was  the  natural  choice  to  utilize  the  MacCormack's 
explicit  finite  difference  scheme^.  For  this  scheme,  the  compu¬ 
tation  for  one  time  step  can  be  accomplished  through  the  following 
simple  functional  relation: 


U 


n+1 
i/ j/k 


L* (At) 


un  .  , 

ii  3/k 


where  the  differencing  operator,  L*(At),  is  given  by 


(41) 


L* (At)  =  L*  (At)  +  L*  (At)  +  L*  (At)  (42) 

£  n  £ 

and 

u"  .  .  =  U  (n*At,i*A£, j*An,k*AC) .  (43) 

This  form  of  unsplit  algorithm  exploits  the  fundamental  character¬ 
istics  of  the  CRAY-1  system,  in  that  it  operates  on  a  single  memory 
path.  This  results  in  a  substantial  saving  on  memory  loading  and 
we  achieve  the  same  final  results  with  further  reduction  in  data 
processing  rate. 
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Since  MacCormack's  scheme  is  conditionally  stable  ,  in  order 
to  avoid  catastrophic  instability,  the  Courant-Fredrich-Lewy  (CFL) 
condition  must  be  satisfied.  For  the  turret  oriented  coordinate 
system,  the  allowable  time  increment  derived  from  the  stability 
analysis  can  be  given  by: 
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where  the  velocity  components  are  defined  as 
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The  CFL  number  is  generally  in  the  raigne  of  0.0  and  1.0.  However, 
a  conservative  choice  of  Courant-Fredrich-Lewy  (CFL =0.8)  condition 
was  implemented  to  have  optimum  time  step  with  stable  numerical 
evolution. 
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We  use  the  solver  to  initialize  the  planes  or  "pages"  in  the 
(n-C)  planes  and  to  sweep,  for  each  time  step,  in  the  C-direction. 
The  vector  length  for  the  present  case  was  62  which  is  also  the 
longest  index  array  of  the  three-dimensional  mesh  system.  The 
computational  facilities  were  utilized  with  remote  access  via 
long  distance  dial-up  and  the  SAMNET  network.  Before  starting 
the  final  computation,  we  solved  the  governing  equations  (6)  through 
(23)  together  with  the  initial  and  boundary  conditions  (27)  through 
(32)  for  the  first  one  hundred  time  steps  at  low  CFL  number.  The 
object  was  to  select  the  proper  mesh  sizes  and  to  estimate  the 
data  processing  rate.  Our  parameters  were  exponentially  stretching 
the  coefficient  along  the  normal  direction  from  the  turret  and 


flat  surfaces  and  the  location  of  the  outer  boundary.  The  results 
are  presented  in  Table  1.  The  best  choice  for  the  coefficients 
was  1.25  (normal  to  the  flat  surface),  5.0  (normal  to  the  turret 
surface)  and  7.0  (the  location  of  the  outer  boundary).  The  finer 
mesh  sizes  were  1.12  *  10  3m  and  0.906  x  I0_3m,  respectively,  and 
the  angular  displacement  between  nodes  was  6°.  These  mesh  sizes 
seem  to  provide  satisfactory  time  step  increment  during  computation. 
Further  refinement  resulted  in  a  stiffen  behavior  of  the  trans¬ 
formation  metric  elements  or  unrealistic  time  step.  For  these 
mesh  sizes,  the  data  prc cessing  rate  improved  as  the  computation 
advanced.  The  final  value  achieved  for  the  data  processing  rate 
was  4.8  *  10  5  seconds. 

Finally,  to  eliminate  the  unwanted  numerical  oscillations  in 
the  computational  domain,  the  artificial  viscosity  terms  are 

included  in  our  solution  procedure.  For  this  purpose,  the  fourth 
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order  pressure  damping  terms  as  rationalized  by  MacCormack 

were  implemented.  In  each  sweep  direction,  these  expressions 

are  of  the  form: 
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The  values  for  the  coefficient  6  were  between  two  and  four.  These 
damping  terms  were  significant  only  in  the  region  of  high  pressure 
oscillations . 

Once  the  initial  phase  of  computation  was  completed  and  some 
anticipated  periodicity  in  the  flow  variables  was  observed,  the 
analysis  of  density  fluctuations  and  wake  structures  was  attempted 
to  compute  the  correlation  coefficient  along  120°  AZA. 
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RESULTS  OF  EXPONENTIALLY  STRETCHING  THE  COEFFICIENT 


SECTION  5 

RESULTS  AND  DISCUSSION 


Before  evaluating  the  various  aspects  of  our  computed 
results,  we  present  a  brief  justification  for  the  first  peak, 
in  the  calculation  of  the  wake  function,  referred  to  in  the 
turbulence  model. 
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The  prime  motivation  behind  improving  Cebeci's  formulation 

was  to  avoid  computing  the  outer  edge  of  the  boundary  layer. 

In  the  present  model'*’*’,  instead  of  solving  the  displacement 

thickness  integral,  the  distribution  of  vorticity  is  used  to 

determine  the  length  scale.  For  the  time  dependent  calculation, 

the  function,  TF  (=  L  •  |  tu  [  )  depends  strongly  on  the  vorticity  which 

sheds  periodically  in  the  downstream  direction.  This  fact  is 

illustrated  in  Figures  6  through  8,  where  the  function  is  plotted 

against  the  nodal  point  location  along  the  normal  direction  of 

the  turret  shoulder.  The  instantaneous-functional  values  at  the 

azimuth  angles  120°  AZA  (Figure  6) ,  180°  AZA  (Figure  7)  and  240° 

AZA  (Figure  8)  are  picked  up  to  demonstrate  the  typical  behavior. 

It  is  certain  that  a  serious  difficulty  is  posed  to  uniquely 

and  consistently  define  TF  (and,  therefore,  L  ) ,  if  the 

rricix  max 

true  maximum  of  TF  is  taken.  As  the  vortices  develop  and  shed, 

the  constant  TFmax  (shown  in  the  figures  by  a  solid  point)  changes 

its  magnitude  and  location  (L  )  abruptly.  This  induces  a  sub- 

ma  x 

stantial  error  in  computing  e  .  In  order  to  avoid  such  irregular¬ 
ity,  we  always  take  first  peak  of  the  function  TF  as  TFmax* 

In  Figure  9  we  given  the  value  of  TF  along  the  azimuth 

angle  0°  through  360°  at  the  turret  shoulder.  It  can  be  observed 

that  in  the  downstream  direction,  for  the  azimuth  angles  90° 

through  270°,  the  function  TF  is  bounded  to  a  reasonable  limit 

max 

even  though  the  time  dependent  variation  exhibits  a  marked  dif¬ 
ference.  (This  limit  was  four-folds  larger  when  the  true  maximum 
was  taken.)  Such  consideration  ascertained  a  convincing  definition 
of  eddy  viscosity  coefficient  c,  as  shown  in  Figure  10.  The 
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The  Function  TF  Along  240°  AZA 


3.0616 


Various  Time  Levels. 
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normal  distance  R  is  normalized  by  the  turret  radius  rt. 
effective  eddy  viscosity  is  within  one  turret  radius  and  remained 
in  this  range  throughout  our  computation. 

During  the  computation  we  recorded  the  pressure  and  Mach 
number  variation  at  twelve  different  locations  in  the  wake  region 
(Figure  11) .  The  points  in  the  plane  were  equally  distributed 
in  the  planes  parallel  to  the  flat  surface  at  heights  Y/R  = 

1.427  x  10  2 (close  to  the  flat  plate)  and  Y/R  =  1.2304  (at  the 
turret  shoulder) .  Figures  12  to  14  give  pressure  variation  at 
some  of  these  locations  for  the  cases  (1)  without  suction;  and 
(2)  with  suction.  It  is  evident  that  the  pressure  oscillations 
reduce  considerably,  in  particular  farther  from  the  turret,  for 
the  flow  with  suction.  This  feature  also  indicates,  as  we  shall 
discuss  later,  that  the  flow  with  suction  should  be  well  organized 
for  the  large-scale  motion. 

For  numerical  solution,  the  comparison  with  experimental 

data  provides  the  validation  of  the  computational  procedure  as 

well  as  guidelines  for  the  subsequent  analysis.  The  experimental 
3  4 

results  '  relevant  to  our  study  are  essentially  parametric  in 
nature  and  offer  the  basic  measurements  at  selected  stations. 

The  pressure  coefficient  distributions  at  the  turret  shoulder 

are  presented  for  the  cases  without  suction  (Figure  15)  and  with 

suction  (Figure  16) .  The  variable  THETA  is  the  azimuth  angle 

measured  i  the  clockwise  direction  from  the  oncoming  free  stream. 

The  different  curves  represent  the  pressure  coefficient  variation 

during  one  principal  time  period.  The  mean  value  agrees  quite 
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well  with  the  experimental  data  .  This  agreement  improves  for  the 
flow  with  suction,  particularly  around  270°  AZA.  It  is  important 
to  note  that  the  presence  of  the  cuto’  •  ,  from  azimuth  angles  96° 
to  144°  has  not  affected  the  pressure  coefficient,  as  it  is  already 
immersed  in  the  separated  flow.  This  M.nts  that  the  cutout  does 
not  influence  the  upstream  separation  mechanism  irrespective  of 
the  flow  conditions  at  the  turret  surface.  In  the  wake  region, 
from  120°  AZA  to  240°  AZA,  the  pressure  recovery  is  faster  and 
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Figure  12.  Pressure  Coefficient  at  -=1.2394 


WITHOUT  SUCTION 


Figure  14.  Pressure  Coefficient  at 


NTAL 


f  pressure  Distribution  for  Flo' 


ssure  Distribution  for  Flow  With  Suction 


the  turret  surface  experiences  pronounced  pressure  variation 
for  the  flow  with  suction.  It  can  be  interpreted  as  if  along 
the  rear  stagnation  line,  there  is  a  tendency  for  the  flow  to 
reattach  and  in  this  vicinity,  the  separation  effects  are  stronger. 

We  also  compare,  in  Figures  17  to  19,  the  Mach  number 
variations  in  the  rear  wake  regions  for  the  flow  without  and  with 
suction.  The  experimental  value  is  available  only  along  the 
azimuth  angle  120°  AZA  at  the  point  A2-  The  characteristic 
time,  TCH,  is  the  ratio  of  the  turret  diameter  and  free  stream 
velocity.  The  mean  local  Mach  number  value  approaches  asymp- 

4 

totically  to  the  time  averaged  experimental  measurements  for 
both  the  cases  (Figure  17).  For  the  flow  with  suction,  the  time 
dependent  variation  as  well  as  the  mean  vali  e  are  higher  than 
their  counterpart.  This  fact  suggests  that  a  small  uniform 
suction  rate  of  0.135  kg/s  eliminates  most  of  the  retarded  fluid 
particles  and  it  intensifies  the  vortical  motion  in  the  wake. 

This  argument  is  further  confirmed  in  Figure  20  which  shows 
the  velocity  field  in  the  turret  shoulder  plane  for  the  flow 
without  suction.  As  the  turret  cutout  is  entirely  within  the 
separated  flow,  it  does  not  cause  any  amplification  to  the  process 
of  vortex  shedding.  On  the  contrary,  a  small  separation  bubble 
appears  to  fill  in  for  the  surface  deflection  and  the  oscillatory 
flow  behaves,  qualitatively,  like  flow  past  a  circular  cylinder2^. 
To  represent  the  wake  structure  in  a  proper  perspective,  a  simple 
interpolation  is  used  to  transform  the  computed  results  from  the 
body  oriented  coordinate  system  (y,8)  to  the  Cartesian  coordinate 
system  (x,z).  To  illustrate  this  procedure,  we  take  four  mesh 
points,  1,  2,  3,  4,  which  are  nearest  to  the  Cartesian  coordinate 
system  point  '5'  (Figure  21).  With  the  help  of  simple  geometry, 
the  curvilinear  coordinates  (y^,0^)  of  the  point  '5'  can  be  found. 
While  the  fluid  properties  (FP)  at  the  points  (y^e^),  (y2,02), 
(73.63)  and  (y 4»®4>  are  known,  by  virtue  of  solving  complete 
Navier-Stokes  equations,  the  fluid  properties  (velocity,  density, 
Mach  number,  etc.)  can  be  computed  at  the  point  (75,65)  by  using 
the  following  relation. 
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Figure  20.  Instantaneous  Velocity  Field  for  Flow  Without 
in  the  Turret  Shoulder  Plane. 
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It  may  be  recalled  that  the  point  '5'  is  a  Cartesian  coordinate 
system  point  and  such  points  can  be  selected  along  the  line 
x  =  constant.  The  important  aspect  remaining  is  to  identify  the 
neighboring  mesh  points  where  the  computation  is  performed.  For 
the  core  without  suction,  the  wake  profiles  and  Mach  number  levels 
are  presented  in  Figures  22  and  23,  respectively. 

When  a  uniform  suction  is  applied,  a  striking  contrast  in  the 
flowfield  is  observed,  as  presented  in  Figures  24  to  26.  This  is 
an  excellent  example  of  the  three-dimensional  flow  interaction  at 
the  turret  shoulder  where  the  inflow  from  the  top  has  sharply 
divided  the  wake  into  two  parts.  The  local  Mach  number  recorded 
in  the  vicinity  of  rear  stagnation  point  was  0.6.  At  this  location, 
the  peridoic  mass  exchange  due  to  unsteady  effect  is  evident  but 
the  characteristic  flowfield  remains  similar.  This  behavior  is 
very  much  in  compliance  with  the  earlier  discussion  regarding  the 
pressure  coefficient  and  Mach  number. 

Complementary  to  the  turret  shoulder  plane  is  the  cross- 
sectional  plane  (0°  -  180°  azimuth  angles)  for  which  the  instan¬ 
taneous  velcocity  field  is  presented  in  Figure  27.  For  the  case 
when  the  turret  is  treated  as  a  solid  shell,  the  upstream  influence 
is  confined  to  the  flat  plate  and  is  within  one  turret  radius. 

In  the  downstream  direction  the  wake  region  is  extended  up  to 
three  times  the  turret  radius  and  is  mainly  affected  by  the 
separated  flow  from  the  turret  top.  The  reverse  flow  and  its 
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Figure  23.  Planar  Plot 


time  dependent  variation  provide  the  mechanism  for  vortex  sheddiny 
and  the  relative  movement  of  reattachment  point.  To  estimate  the 
three-dimensional  Mach  number  distribution  close  to  the  hemi¬ 
spherical  turret  surface,  we  make  use  of  the  medium  plane  anyle 
(MPA).  The  0°  MPA  corresponds  to  the  turret  shoulder  while  90°  MPA 
corresponds  to  the  turret  top.  Figure  28  gives  one  such  instan¬ 
taneous  distribution  for  the  same  time  level  as  that  of  the  vel¬ 
ocity  field  shown  in  Figure  27.  It  is  interesting  to  observe  that 
the  maximum  Mach  number  occurs  at  89°  AZA  and  276°  AZA  and  the 
minimum  at  the  rear  stagnation  line  neighborhood.  This  phenomena 
demonstrates  a  strong  reverse  flow  on  the  turret  surface  and  that 
the  flow  diverges  from  the  stagnation  line.  A  marked  difference 
in  these  characteristics  is  observed  in  Figures  29  and  30  when 
the  turret  is  considered  to  a  porous  shell.  As  a  result  of  the 
small  uniform  suction,  the  separated  and  reverse  flow  regions  are 
reduced  substantially.  In  this  case  the  flow  at  the  turret  top 
is  attached  and  the  maximum  Mach  number  is  located  at  192°  AZA. 

The  local  Mach  number  close  to  the  surface  crosses  the  sonic  limit 
after  42°  MPA. 

A  significant  insight  into  the  three-dimensional  separated 

flow  is  offered  by  the  limiting  streamline  pattern,  which  are, 

in  the  limiting  case,  the  skin  friction  lines.  These  lines 

satisfy  the  topological  criteria  for  the  singular  points.  This 

confirmation  was  provided  in  our  earlier  work^  for  the  turret  as 

well  as  for  the  flat  plate  surfaces.  For  this  study  also,  the 

primary  separation  line  and  the  horseshoe  vortex  structure  for 

large  protuberance  (Figures  31  and  32)  show  an  excellent  agreement 

9 

with  the  work  of  Sedney  and  Kitchens  .  In  Figures  33  and  34,  we 
present  two  typical  limiting  streamline  patterns  (for  one  principal 
time  period)  on  the  turret  surface.  It  can  be  observed  that  for 
the  flow  without  suction,  the  turret  top  experiences  a  highly 
random  flowfield  behavior  and  the  turbulence  is  prominent  in  the 
cutout  side.  This  "opening"  switches  from  cutout  side  to  smooth 
side  and  back  as  the  computation  progresses.  Thus,  the  separation 
line  encounters  an  oscillatory  phenomena  beyond  60°  MPA.  This 
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Figure  33.  Instantaneous  Limiting  Streamline  Pattern  on  the 
Turret  Surface  for  Flow  Without  Suction. 
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undesirable  process  is  completely  eliminated  for  the  flow  with 
uniform  suction  where  the  limiting  streamline  pattern  exhibits 
a  well  organized  motion.  The  disturbances  caused  by  the  turret 
cutout  are  only  local  and  they  do  not  affect  the  separation  line 
location,  in  general. 

In  order  to  describe  the  oscillatory  flowfield  structure 
along  120°  AZA  at  the  turret  shoulder,  the  statistical  properties 
of  unsteady  density  are  analyzed.  For  compressible  turbulent 
flow,  the  time  mean  and  root-mean-square  (RMS)  values  were  cal¬ 
culated  as 


P  = 


If 
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where  T  =  0.0056  second  is  equal  to  two  principal  time  periods 
P 

and  is  a  function  of  Strouhal  number,  turret  diameter  and  the  free 
stream  velocity.  Since  the  time  increment  for  the  present  analysis 
is  very  small,  the  straightforward  integral  procedure  (48)  and  (49) 
were  adopted.  The  density  values,  p_.  at  grid  point  location  j, 
were  recorded  for  3400  consecutive  time  iterations  (CFL  =  0.8)  and 
were  integrated  for  both  cases.  The  results  are  shown  in  terms 
of  mean  density,  Figure  35,  and  <p'>/<p'>  ,%,  Figure  36,  as 

functions  of  distance  AR/R  from  the  turret.  The  mean  density  near 
the  turret  surface  reduces  when  suction  is  applied,  indicating 
rapid  expansion,  but  recovers  its  free  stream  value  within  one 
turret  radius.  It  was  also  observed  that  during  the  time  period, 

T  ,  the  difference  in  the  mean  fluctuation  values  of  density  be¬ 
tween  the  flow  without  and  with  suction  case  were  negligible. 
However,  the  RMS  density  fluctuations  are  higher  for  flow  with 
suction.  This  was  true  for  experimental  as  well  as  computational 
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Figure  36.  Comparison  of  RMS  Density  Fluctuations  with 
Experimental  Data. 


values  predicting  a  thinner  shear  layer.  By  virtue  of  clustering 
the  mesh  points  near  the  turret  surface,  there  seems  to  be  a 
precise  computational  information  in  terms  of  density  data.  In 
this  region,  wo  did  not  have  any  experimental  value  to  compare 
with.  Away  from  the  turret,  even  though  the  functional  variation 
is  similar,  the  computed  results  underpredict  the  density  shear 
layer.  For  the  experimental  data,  the  ratio  of  the  two  peak  RMS 
density  values  was  3.6780,  while  for  computed  results  it  is  2.850. 
This  discrepancy  can  be  partially  attributed  to  the  difference  in 
the  test  and  simulation  conditions  and  the  possible  excessive 
vibration  of  the  model  during  the  test  procedure^.  The  computed 
results  also  can  be  improved  by  reducing  grid  space  and  intro¬ 
ducing  more  mesh  points.  This  might  tax  severely  the  available 
CRAY-1  memory  and  the  computer  resources.  It  was  not  possible  at 
the  present  time,  but  a  rigorous  study  in  the  three-dimensional 
mesh  generation  area  is  one  of  our  future  goals. 


Continuing  our  analysis,  the  correlation  coefficient, 
for  the  density  data  at  mesh  point  location  j  and  j  was 

1  c  Jm  Jn 

evaluated  as 


R  , 
mn 


(50) 


m,n  =  1,2, •••30 


Evidently,  R  is  bounded  and  R  =1.0  indicating  that  the  density 
1  mn  mm  J 

values  at  the  mesh  point  location  j  are  correlated.  We  have 

m 

plotted,  in  Figures  37  and  38,  the  correlation  coefficients  for 

the  first  six  mesh  locations  to  emphasize  the  near  field  structure. 

During  the  time  T  ,  for  the  flow  without  suction,  the  function 

Rmn  does  not  show  any  correlation  and  the  shear  layer  is  highly 

unstable.  The  uniform  suction  provides  a  coherence  in  the  fluid 

flow  behavior  and  a  specific  definition  of  shear  layer  width. 

Beyond  the  shear  layer,  the  numerical  value  of  R  oscillates  close 

mn 

to  zero. 
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Finally,  along  120°  AZA,  a  spectral  analysir  is  done  for  all 

possible  modes  of  oscillations.  Such  evaluation  is  one  of  the 

important  aspects  of  unsteady  fluid  dynamics.  We  accomplish  this 

goal  through  the  ninety  term  Fourier  series.  The  computation  is 

done  for  all  the  grid  points  and  Figure  39  is  the  representative 

graph  for  the  location  AR/R  =  0.486.  The  first  peak  is  observed 

at  the  frequency  349.83  Hz  which  corresponds  to  the  Strouhal 
25 

number  0.2265.  The  higher  modes  of  the  spectral  analysis  exhibit 
a  complete  random  behavior;  no  single  discrete  frequency  is 
dominated  over  the  other.  Based  on  the  phase  angle  analysis 
(Figure  40) ,  we  conclude  that  the  highest  credible  frequency  range 
of  our  numerical  analysis  is  around  6.4  KHz.  This  limitation  can 
only  be  eliminated  by  using  a  much  finer  mesh  point  system. 


AMPLITUDE  (LOG) 


Figure  39.  Spectral  Analysis  Along  120°  AZA. 
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SECTION  6 
CONCLUSIONS 

For  a  Mach  number  of  0.6  and  Reynolds  number  of  12.8  *  106/m, 
the  numerical  simulation  of  time  dependent  three-dimensional 
separated  flow  is  accomplished.  It  is  established  that  a  uniform 
suction  at  the  turret  surface  causes  a  significant  variation  in 
the  wake  structure.  Qualitative  agreement  is  obtained  for  the 
general  flowfield  including  RMS  density  fluctuation  and  correlation 
coefficient.  However,  in  order  to  achieve  better  resolution,  the 
optimum  rearrangement  of  grid  points  in  all  three  directions 
should  continue  to  be  stressed  in  future  investigations. 
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APPENDIX  -  VORTICITY  FUNCTION 


The  vorticity  in  the  three-dimensional  space  is  given  by 

co  =  V  x  u 
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where  fj^,  fi2,  are  the  unit  vectors  in  x,  y,  and  z  directions, 
respectively.  The  magnitude  of  this  vorticity  vector  is,  then, 
given  by 
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To  evaluate  the  magnitude  of  vorticity  in  the  body  oriented 
coordinate  system  5/  n»  C,  we  use  the  Jacobian  of  transformation, 
Jn,  to  give 
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On  the  turret  surface 
^-direction  vanish.  Thus, 
expressed  as 


the  derivatives  in  the  ^-direction  and 
the  surface  vorticity,  ww,  can  be 
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